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Higher-dimensional models of neutrino physics with one or more right-handed neutrinos in the bulk 
have attracted considerable attention in recent years. However, a critical issue for such models is to 
find a way of introducing the required flavor dependence needed for generating neutrino oscillations. 
In this paper, we point out that a natural "minimal" framework that accomplishes this can be 
constructed by combining the bulk-neutrino hypothesis for right-handed neutrinos with the split- 
fermion scenario for left-handed neutrinos. This combination leads to a unique flavor signature for 
neutrino phenomenology which easily incorporates large flavor mixing angles. This hybrid scenario 
also has a number of additional important features. For example, one previous difficulty of the split- 
fermion scenario applied to neutrinos has been that the mass matrix is exponentially sensitive to 
neutrino displacements within the brane. However, in our hybrid scenario, the interactions between 
the brane and bulk naturally convert this dependence from exponential to linear. Another important 
feature is that our hybrid scenario provides its own natural regulator for Kaluza-Klein sums. Thus, 
in our scenario, all Kaluza-Klein summations are manifestly finite, even in cases with multiple extra 
dimensions. But most importantly, our mechanism completely decouples the effective neutrino flavor 
mixing angles from the sizes of the overlaps between the neutrino wavefunctions within the brane. 
Thus, we are able to obtain large neutrino mixing angles even when these neutrinos have significant 
spatial separations and their overlaps vanish. 



I. INTRODUCTION AND OVERVIEW 

Many physicists consider neutrino oscillations to be 
the first clear signature for physics beyond the Standard 
Model. With the most recent data from SNO and Super- 
Kamiokande, there remains virtually no doubt about the 
existence of neutrino flavor oscillations. Indeed, the pa- 
rameter space for neutrino mass differences and mixing 
angles has already become significantly constrained. 

There remains, however, the paramount issue con- 
cerning how these experimentally determined parame- 
ters can be accommodated within a theoretically moti- 
vated model. While there have been numerous avenues 
that have been explored within recent years, models with 
large extra dimensions have received considerable atten- 
tion ever since it was realized that extra 
spacetime dimensions have the potential to provide new, 
intrinsically geometric perspectives on the hierarchies be- 
tween various energy scales present within and beyond 
the Standard Model (SMV These include the Planck 
scale , the GUT scale , and the string scale , 
as well as a possible scale for SUSY-breaking 0. Neu- 
trino phenomenology has also been investigated within 
this higher-dimensional context, with the fundamental 
idea 0,13 being that since the right-handed neutrino is a 
Standard-Model singlet, it need not necessarily be bound 
to the brane to which the other Standard-Model particles 
are restricted. The right-handed neutrino can therefore 
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propagate into the bulk of extra spacetime dimensions 
and accrue an infinite tower of Kaluza-Klein (KK) exci- 
tations, all of which will then take part in neutrino oscilla- 
tions [3- This idea, along with numerous variations, has 
spawned a relatively large literature investigating vari- 
ous aspects of higher-dimensional neutrino phenomenol- 
ogy (see, e.g., Refs. 00,10). 

The challenge for these bulk-neutrino models, however, 
is to accommodate a suitable flavor structure which can 
trigger the observed neutrino oscillations. To date, two 
different directions have been explored in the literature. 

One direction involves generalizing the original idea of 
Refs. by introducing a separate right-handed bulk 
neutrino for each flavor of neutrino on the brane. This 
process thereby extends flavor into the bulk. However, 
in such cases, the brane/bulk couplings become arbitrary 
3x3 mixing matrices whose parameters are in principle 
undetermined. Moreover, the three bulk neutrinos can in 
principle correspond to different extra spacetime dimen- 
sions with different radii. Thus, one obtains a scenario 
with many undetermined parameters governing neutrino 
masses and mixing angles. Indeed, to a large extent, 
the role of the extra dimensions in such scenarios is re- 
duced to providing a set of Kaluza-Klein states which 
function (from a four-dimensional perspective) as addi- 
tional sterile neutrinos. In contrast, flavor neutrino os- 
cillations continue to be triggered through mixing angles 
which are ultimately introduced by hand, just as in four 
dimensions. This therefore fails to yield any further in- 
sight regarding the structure or origins of the neutrino 
mixings. 

A second, more "minimal" idea first advanced in 
Ref. is to introduce only one bulk neutrino, and to pro- 
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vide this single bulk neutrino with flavor-universal cou- 
plings to all three brane neutrinos. Thus, in so doing, one 
is essentially considering flavor to be a feature internal to 
the Standard Model, a feature which is restricted to the 
brane and which therefore does not extend into the bulk. 
It is even possible to take the brane theory to be flavor- 
diagonal, thereby avoiding mixing angles completely; in- 
deed, in the model of Ref. [9j, the only flavor-sensitive 
feature that distinguishes these brane neutrinos is their 
differing bare Majorana masses m^. Nevertheless, one 
finds 0] that sizable neutrino flavor oscillations between 
the different neutrinos on the brane arise as a result of 
their indirect mixings with the Kaluza-Klcin modes of 
the bulk right-handed neutrino. Thus, in such a model, 
the flavor oscillations are entirely "bulk-mediated" and 
no mixing angles are needed at all: it is the presence of 
the higher-dimensional bulk which is completely respon- 
sible for inducing the flavor oscillations on the brane. 

Although it is possible to consider cases in which the 
couplings between the brane neutrinos and the bulk neu- 
trino are large, it is convenient (and perhaps also phe- 
nomenologically necessary) to consider a so-called "per- 
turbative" limit in which these couplings are relatively 
small. In such cases, it is possible to integrate out the 
effects of the Kaluza-Klein bulk neutrinos in order to ob- 
tain an effective mixing matrix involving only the light 
modes in the theory. Typically, these light modes consist 
of the three brane neutrinos as well as the zero mode(s) 
of the bulk neutrino(s). Because this effective mix- 
ing matrix completely encapsulates the resulting flavor- 
dependence of the model, it inevitably lies at the center 
of any comparison between a given theoretical neutrino 
model and experimental data. Indeed, experimentalists 
will eventually provide a unique numerical mixing matrix 
which incorporates the entire observed phenomenology of 
neutrino oscillations. It will then be the job of the theo- 
rist to solve the "inverse" problem of deducing the set of 
viable underlying neutrino models which can lead to this 
matrix. 

The non-minimal models with many bulk neutrinos 
contain many mixing angles and can therefore give rise 
to relatively diverse effective mixing matrices. However, 
one difficulty with the minimal model of Ref. is that 
it leads to a rather rigid form for this effective mixing 
matrix. Since the minimal model involves only one bulk 
neutrino, it yields an effective 4x4 mixing matrix which, 
as we shall see, has a texture of the form 
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where the first three rows/columns correspond to the 
brane neutrinos and the final row/column corresponds 
to the bulk zero mode. Here X and m are parameters 
associated with the bulk of the higher-dimensional theory 
and its coupling to the brane, while mi are the bare Majo- 
rana masses of the neutrinos on the brane. While we see 



that the bulk physics (through the term X) is responsible 
for yielding a non-diagonal 3x3 brane mixing submatrix, 
this matrix is relatively rigid, with all off-diagonal entries 
forced to be exactly equal. Indeed, a mixing texture of 
this form may have difficulty accommodating observed 
neutrino oscillations of the correct sizes and magnitudes. 

It is the purpose of this paper to propose an alterna- 
tive "minimal" model which leads to a richer effective 
mixing matrix without introducing additional parame- 
ters or sacrificing any of the minimality of the single 
bulk-neutrino scheme. Rather than incorporate a flavor 
structure through differing bare Majorana masses Wj on 
the brane as in Ref. we shall instead incorporate a 
non-trivial flavor structure within the context of the so- 
called "split-fermion" scenario [ljj. As i s weu known, the 
split-fermion scenario represents an intrinsically higher- 
dimensional method of explaining flavor hierarchies for 
the charged matter content of the Standard Model. How- 
ever, as we shall discuss below, this method faces cer- 
tain unique difficulties when attempting to address the 
flavor structure of the neutrino sector. Thus, by com- 
bining the bulk-neutrino scenario with the split-fermion 
scenario, we are able to obtain a "hybrid" minimal model 
which has both strong theoretical motivations as well as 
rich prospects for neutrino phenomenology. 



II. SPLIT FERMIONS 

In the split-fermion scenario |10| ] , the fermions of the 
Standard Model are located within a "fat" brane but cen- 
tered around different positions within the brane. Their 
spatial extent within the fat brane is modeled by a Gaus- 
sian distribution with a typical width er, where a is ap- 
proximately one order of magnitude smaller than the 
width of the brane. The gauge fields and the Higgs are 
assumed to be equally distributed over the width of the 
brane. To obtain the effective four-dimensional couplings 
between different particles, the extra dimensions have to 
be integrated out. Since the overlap of the Gaussian 
wavefunctions of the particles can be very small, the re- 
sulting effective couplings can be extremely suppressed. 
Indeed, it is easy to obtain a suppression by thirty orders 
of magnitude by delocalizing particles at locations which 
are separated by about ten times their Gaussian widths. 

This setup yields a solution to several problems. It can 
be used to suppress the proton decay, to suppress flavor- 
changing operators, or to explain the observed mass hier- 
archy in the lepton sector. Many possible configurations 
of particles inside the brane are possible; some of these 
are discussed in Ref. 0, 0, ^| . 

Relative shifts between left-handed Standard-Model 
doublets and right-handed singlets have been successfully 
used to address the mass hierarchy of quarks and charged 
leptons. But for the same reason the split-fermion sce- 
nario is easy to implement for charged fermions, it leads 
to difficulties when one tries to incorporate neutrinos |14| . 
Experiments suggest with an increasing degree of cer- 
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tainty that neutrino mixing is maximal in the case of 
v T <-> j/„ and almost maximal in the case of v e <-> v T . 
The common textures of the zeroth-order neutrino mass 
matrix that are compatible with such almost-bimaximal 
mixing can be classified by the hierarchy type into differ- 
ent categories — normal, inverted, or degenerate. How- 
ever, all of these situations require several entries to be 
of the same order of magnitude. 

Unfortunately in the split-fcrmion scenario, the cou- 
pling between any two particles is extremely sensitive to 
their relative distance within the brane. The observed 
neutrino mixing can therefore be achieved only by care- 
fully choosing the central locations of the neutrinos inside 
the brane. In addition, suitable positions for the charged 
leptons must be found in order to avoid excessive flavor- 
changing processes. Together, these requirements pu t 
severe constraints on the allowed fermion locations |14| . 
While all experimental constraints can ultimately be ac- 
commodated, the resulting particle map is extremely sen- 
sitive to small perturbations. Thus, even though these 
models reduce the amount of fine-tuning for the Yukawa 
couplings, a considerable degree of fine-tuning continues 
to be necessary. 

There are, of course, various options for ameliorating 
this situation. For example, introducing Majorana neu- 
trinos rather than Dirac neutrinos has been proposed 
as a method of achieving reasonable neutrino masses 
within the split-fermion scenario. The neutrino sector 
of the s plit -fermion scenario has also recently been ex- 
amined llg within the context of the Randall- Sundrum 
model 

In this paper, we will consider a different approach to- 
wards the neutrino sector by asking what consequences 
arise if the two scenarios — split fermions and bulk neu- 
trinos — are examined in a joined framework. In this 
way, we shall be combining the strengths of each indi- 
vidual scenario: the split-fermion mechanism will trigger 
the flavor structure of the neutrino sector, while the bulk 
sector will ameliorate the fine-tuning issues and naturally 
lead to large mixing angles. Indeed, it was already noted 
in Ref . that the combination of having brane neutrinos 
at different locations in the extra dimension would trigger 
a higher-dimensional seesaw mechanism, although this 
idea was not pursued and no consistent framework was 
provided. In this paper, we shall see that split fermions 
provide a natural context for this phenomenon, leading 
to flavor-dependent mixing properties which can be used 
to distinguish between the brane neutrinos. 

There are also other benefits to combining these sce- 
narios. For example, as we have mentioned, one diffi- 
culty of the split-fermion scenario applied to neutrinos 
has been that the mass matrix is exponentially sensitive 
to neutrino displacement within the brane. However, in 
our hybrid scenario, we shall see that the interactions 
between the brane and bulk naturally convert this de- 
pendence from exponential to linear. Thus, the previous 
exponential sensitivity is entirely eliminated. Another 
important feature is that our hybrid scenario provides its 



own natural regulator for Kaluza-Klein sums. Thus, in 
our scenario, all such Kaluza-Klein summations are mani- 
festly finite, even in cases with multiple extra dimensions. 
As we shall see, this arises because the heavy Kaluza- 
Klein modes with wavelengths that are small compared 
with the widths of the brane Gaussians will average out 
when folded with the Gaussian distribution. This thereby 
eliminates the necessity of introducing a cutoff by hand. 

But most importantly, we shall find that our mech- 
anism completely decouples the effective neutrino fla- 
vor mixing angles from the magnitudes of the neutrino 
overlaps within the brane. Thus, we are able to ob- 
tain large neutrino mixing angles between neutrinos, even 
when these neutrinos have rather large spatial separa- 
tions within the brane. 



III. THE HYBRID FRAMEWORK 

We now describe the framework for our hybrid model. 
Clearly, parts of this model will be similar to the model 
considered in Ref. [9|. For concreteness, we will work in 
this section with one bulk neutrino and one compactified 
extra dimension of radius R, although the model can be 
trivially extended to more bulk fields and to more addi- 
tional dimensions. 

The primary features of this model can be summarized 
as follows. We imagine a single extra dimension of radius 
R, as well as a brane of width R' <C R. 

• On the brane, we introduce nf left-handed neutri- 
nos v a \ these are our flavor eigenstates. We as- 
sume that these left-handed neutrinos have no cor- 
responding Majorana masses and are restricted to 
lie within the brane. 

• In the bulk, we introduce a single Dirac fermion 
\& which does not carry any flavor indices and is 
therefore completely flavor-neutral. 

• For the sake of minimality, we assume that the 
left-handed neutrinos do not mix or couple to each 
other directly. Thus, all mixing angles on the brane 
are set to zero. 

• We introduce a single flavor-blind brane/bulk cou- 
pling M* between each of the brane neutrinos v a 
and the bulk neutrino \P. 

• Finally, each of the active brane neutrinos v a is 
assumed to be centered around a different trans- 
verse location within the brane. It is this feature, 
and only this feature, which is introduces a flavor- 
dependence into our model. 

We shall take our spacetime coordinates as x A = 
(x^, y) where are the four uncompactified coordinates 
on the brane and y is the coordinate of the fifth dimen- 
sion perpendicular to the brane. Given the above as- 
sumptions, the action for our model then takes the form 
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S = S v + Sb + S c + h.c. where S v describes the physics 
on the brane, Sh describes the physics in the bulk, and 
•Sc describes the brane/bulk couplings. These individual 
contributions take the form 



S v = 
Sb = 

S c = 



d 4 xdy y~] u a (x, y^a^d^Ug (x, y) 

a=l 

J d A xdy^{x,y) iY A d A ^(x,y) 
J d^xdy ^2 v a (x, y) [ip+(x, y) + ^-(x, y)] 



a=l 
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where Y represent the five-dimensional Dirac matrices. 
Here the index a = 1, .., ri{ runs over all rif flavors on the 
brane. 

In the above, we have chosen to work in the Weyl basis 
in which the Dirac fcrmion ^> can be decomposed into two 
two-component spinors: $ = V'-) T - The bar on 
tp— indicates that ip- transforms as a different Lorentz 
representation than ip+. We refer to ip+ and %j>- as left- 
handed and right-handed components, respectively. We 
have also written tr M = (1, — <r l ), <r M = (1, +cr l ) where a 1 
are the Pauli matrices. 

Given the action in Eq. (J2J), the next step is to com- 
pactify the fifth dimension in order to obtain an ef- 
fective four-dimensional theory. To do this, we shall 
make the following assumptions. First, we shall iden- 
tify y ~ y + 2ttR and impose the orbifold relations 
i>+{-y) = ip+{y) and tf-i-v) = -tf-iv), where y_is 
the coordinate of the fifth dimension. Thus, ip + and ip- 
have the Kaluza-Klein mode expansions 
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Note that as a result of our orbifold boundary conditions, 
does not have a zero mode. We shall also take our 
brane neutrinos to have wavefunctions of the form 



v a {x,y) = G(y - y a ,a)v a (x) 



(4) 



where v a {x) denotes the usual four-dimensional part of 
the fermion wavefunction on our brane and where G(y, a) 
denotes a normalized Gaussian wavefunction centered 
around location y = with width a: 
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Note that for simplicity, we are assuming that each brane 
neutrino has the same width a <C R- 



Inserting the wavefunctions in Eqs. 10 and J2J into 
the action J2J) and integrating over the fifth dimension, 
we obtain the effective four-dimensional actions: 
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Note that the fields are now functions of x M only, 
and m and m™ ± are volume-suppressed n-dependent 
brane/bulk couplings resulting from the rescaling of the 

individual ip^L , ip^ Kaluza-Klein modes: 
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As we see, the brane/bulk coupling is suppressed by a 
volume factor \J a/ R. Moreover, for non-shifted brane 
neutrinos with y a = 0, we see that the couplings m" _ 
actually vanish. Thus, it is only the possible displace- 
ment of the antisymmetric fermion modes 'ip^f 1 away from 
y a = which permits them to couple to the bulk zero 
mode. 

Given the Lagrangian in Eq. @, we immediately see 
that the SM flavor-eigenstate neutrinos v a on the brane 
will mix with the entire tower of Kaluza-Klein states of 
the higher-dimensional ^ field, even though they do not 
mix directly with each other. Defining 
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we see that the mass term of the Lagrangian © takes 
the form Af T A4AT where the mass matrix M. is given by 
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Of course, in this notation, the first row/column repre- 
sents Hf independent rows/columns, one for each flavor 
neutrino on the brane. The fact that this mass matrix 
is non-diagonal implies that each of the rif brane neutri- 
nos will undergo oscillations with the bulk KK modes. 
Through these mixings with the bulk KK neutrinos, the 
brane neutrinos will thereby undergo effective flavor os- 
cillations with each other. 

The next step is to determine the eigenvalues and 
eigenvectors of the mass matrix ©. We can begin 
this task as follows. Defining the eigenvector v = 
(v a ,v(°\v+\v_ , v+\ ...), we see that the eigen- 
value equation Mv — Xv yields the simultaneous equa- 
tions: 
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Combining the last two of these equations yields a partial 
solution for u± , 



Note that for numerical purposes, Eq. i|13|l provides a 
particularly useful simplification of the original infinite- 
dimensional eigenvalue problem. Solving Eq. I|13() for A 
and v and inserting these results into Eq. then en- 
ables one to obtain the full infinite-dimensional eigenvec- 
tor v. 

Thus far, we have not made any approximations. How- 
ever, in order to proceed analytically, we shall now take 
the so-called "perturbative" limit mR <C 1 in which the 
brane/bulk coupling is small compared with the com- 
pactification scale. Introducing this limit is one way of 
ensuring that the brane sector of our theory (or, as we 
shall see, a subset of that sector) will experience little or 
no net loss of probability into the "sterile" tower of bulk 
excitations. In this limit (see, e.g., Ref. Q), the eigen- 
values of the excited KK neutrino modes will be approx- 
imately ±n/R, corresponding to the bulk/bulk sector of 
the mixing matrix, while the remaining rif + 1 eigenval- 
ues will be much smaller than the masses of the excited 
modes. Specifically, we see from the form of the reduced 
mass matrix in Eq. (|14|l that the eigenvalues of these 
remaining light modes will scale either as A ~ M 2 a or 
as A ~ rn. However, the assumption that mR <C 1 im- 
plies that XR <C 1 in either case: if A ~ m then clearly 
XR ~ mR < 1, while if A ~ M 2 a then 



XR ~ M 2 aR ~ (mi?) 2 < 1 
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and inserting this into the top line of Eq. l|10fl along with 
the definitions in Eq. J2J) yields the relation 
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Thus, we see that the task of solving the remaining 
equations has now been reformulated as an equivalent 
(rif + 1) x (rif + l)-dimensional eigenvalue problem of the 
form 



M(X)v = Xv 



(13) 



where v = (v a ,v^) and where Ai(X) is a A-dependent 
"reduced" mass matrix of the form 
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Here m a p = (M 2 a/ir)f a {3, where the dimensionless coef- 
ficients fa/3 (A) are given by 
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for these remaining light modes. As wc shall demonstrate 
in Appendix A, this result may also be derived through 
a direct perturbative diagonalization of our original mass 
matrix. 

Note that the only case in which it is not justified to 
pass from Eq. l(TS|l to Eq. (|T7|) occurs when y 2 is exactly 
equal to — yi, for in this case the second term in Eq. H15|) 
vanishes exactly for all n and is thus always smaller than 
the first. However, since A is significantly suppressed for 
the light modes we are considering, we see that the re- 
sulting value of f a p in this case is significantly smaller 
than it is for other values of j/i and yi. The approxi- 
mation in Eq. I|17|) . which sets this value of f ol p to zero, 
therefore continues to be roughly accurate even in this 
case, indicating a strongly suppressed value for f ap in 
the region near yi — —y\. 

It is relatively straightforward to evaluate the KK sum 
in Eq. (|T7|) , especially since the presence of the Gaussian 
width for the brane neutrinos serves as a natural regu- 
lator for the KK sum which renders it manifestly finite. 
Introducing the variable k = (a/R)n and noting that 
a <C R, we see that the sum can be approximated as an 
integral, yielding 

f°° dk 

fap ~ - J cxp(-fc 2 /7r) sin[fc(y Q +yp)/<j] 
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Note that the error function is an odd function, 
Erf (— x) = —Erf (a;), with approximate behavior 



Erf (x) 



x for |x| < 1 

sign(x) for |x| > 1 
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Thus, defining the rescaled (dimensionless) displace- 
ments y a = y/ny a /(2a), we conclude that for the light 
modes in our theory, our reduced (rif + 1) x (rif + 1)- 
dimensional mass matrix in Eq. I|14|l takes the form 
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•27a) for |27 a +27/3I £ 1 • 
(21) 

We see from Eq. (|2*TJl that when the neutrino displace- 
ments are of the same order as their Gaussian widths, 
the entries of this mass matrix depend linearly rather 
than exponentially on the displacements! This is an im- 
portant observation, implying that the usual exponential 
sensitivity to the fermion displacements has been elimi- 
nated in our hybrid scenario, with the interaction with 
the bulk neutrino serving to convert this sensitivity from 
exponential to linear. Moreover, we also observe the in- 
teresting fact that the interactions with the bulk neutrino 
actually remove all sensitivity to the neutrino displace- 
ments when these displacements become significantly 
larger than their Gaussian widths, with mass-matrix val- 
ues saturating near the universal value m a p m —~M%cr 
in such cases. Thus, in such cases, all flavor dependence 
is lost. 

We also notice from Eq. I|21|l that our effective mixings 
on the brane are not translationally invariant within the 
brane, since they depend not on the difference y a — yp, 
but rather on the sum y a + yp . Indeed, our off-diagonal 
mixing terms cancel only for y a = —yp rather than 
Va — yp- Moreover, since this dependence is on the sum 
rather than the difference, the net distance from the ori- 
gin at y — becomes an important quantity, determining 
when the error function passes from its linear to constant 
regime. Once again, however, these features represent the 
effects of the mixings with the bulk neutrino modes, for it 
is the KK orbifold relations just above Eq. 10 which im- 
plicitly select the origin y = as a special point, thereby 
breaking translational invariance within the brane. 

These features imply that it is possible to obtain large 
effective flavor mixing angles for the active neutrinos on 
the brane even when the overlaps between these neutrinos 
vanish. In other words, large neutrino overlaps are no 
longer needed for large effective mixing angles. 

To illustrate this behavior more concretely, let us con- 
sider a simplified two-flavor example in which we have 



only one effective mixing angle. Given the result in 
Eq. i|21|) . we may identify this effective mixing angle 8 
through the relation 



tan26> 
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Erf(27 2 )-Erf(2y 1 
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For a given fixed value of 2/1 (the position of the first 
neutrino), we can now examine how this effective mixing 
angle depends on 2/2 (the position of the second neutrino). 

There are several cases to consider. If |7i | <§C 1 (so that 
the first neutrino is centered very close to the origin), we 
find that tan 28 — > 2 as 2/2 - > ±00. Thus, we find that 
sin 2 28 — > 4/5 as 2/2 — > ±00, which represents nearly max- 
imal mixing. Note that this result applies even though 
the second neutrino is being pulled relatively far from the 
first. Indeed, one can verify the mixing remains nearly 
maximal in this case regardless of the value of 2/2 • 

By contrast, if \yi\ > 1, our asymptotic behavior 
changes. For concreteness, let us take yi positive (oth- 
erwise we can change all y — > —y without altering the 
results). In this case, Erf (2/1) 1 and we see from 
Eq. 1)22(1 that tan 28 — > 00 as 2/2 — * 00 while tan 28 — + 2 as 
2/2 — * —00. Thus, asymptotically, we have relatively large 
mixing angles in either case. However, near 2/2 ~ — 2/1 j 
we find that our mixing angle passes through zero and 
changes sign. Thus, there is a region near 2/2 ~ —yi 
for which smaller mixing angles can be realized without 
exponential fine-tuning. 

Finally, for |7i| 1, the resulting behavior interpolates 
between the two behaviors described above. Again taking 
2/1 > for simplicity, we find that sin 2 28 asymptotes to 
values between 1/2 and 4/5 as 2/2 — * —00, and between 
4/5 and 1 as 2/2 — + +00. However, as in all cases, the 
effective mixing angle hits zero at 2/2 = — 2/1 ■ 

This behavior is shown in Fig. ^ for 2/1 = 0, 1, 2, 3. In 
each case, we observe that large mixing angles dominate, 
arising in most regions of the parameter space. Never- 
theless, we see that small and intermediate mixing angles 
are easy to achieve without significant fine-tuning in the 
region near 2/2 = —yi- 

As indicated in the Introduction, the flavor mixing in 
this model is significantly different from that of an earlier 
model in which all brane neutrinos are located at y = 
with zero width, distinguished only through their differ- 
ent "bare" Majorana masses m.j on the brane. In such a 
model, the brane neutrinos couple to only the even bulk 
modes t/u!^ , whereupon we see (following the same analy- 
sis as performed here) that the effective (rif + 1) x (rif + 1) 
"reduced" mixing matrix A^(A) takes the form shown in 
Eq. where 



X 



A 



7rm 2 i?cot(7ri?A) 



(23) 



The important point here is that X is completely flavor- 
independent. Thus all of the off-diagonal entries in the 
effective mass matrix are identical, leading to a relatively 
rigid flavor structure. 
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FIG. 1: Effective flavor mixing angle sin 2 26 as a function of 
neutrino positions ( J/1,2/2) in the case of two active flavors on 
the brane. We plot sin 2 26 as a function of j/2 for yi =0,1,2,3 
[curves (a) through (d) respectively]. In each case, we see 
that large mixing angles are fairly common, with small mixing 
angles emerging only near j/2 ~ — yi. 



Thus far, we have focused on the case with only a single 
extra spacetime dimension. It is therefore worth under- 
standing how the relevant mass scales in the mass ma- 
trix i|2U|) are altered if we consider the analogous situation 
with S > 1 extra dimensions. With 6 extra dimensions, 
our three-brane has a "width" in 6 directions, so that 
our brane neutrinos are 5-dimensional Gaussians. Trac- 
ing our previous steps, this implies that the brane/bulk 
coupling m defined in Eq. (JJJ will accrue a further volume 
factor of yj a / R for each additional spacetime dimension, 
so that we can generally replace 



/?? 



mia/R^ 5 - 1 ^ 2 ~ M*(a/R) s/2 . (24) 



On the other hand, it is shown in Appendix B that the 
entries of m a/ 3 continue to scale as M 2 a regardless of the 
number or radii of extra dimensions. Thus, for general 
S, we obtain an effective (nt + 1) x (rif + l)-dimensional 
mass matrix whose entries scale as 



M 



M*(7 f al 3 

{a/R) 5 ' 2 



(a/R) 5 / 2 




(25) 



where f a p denotes a ^-dependent set of flavor-dependent 
coefficients of order one. We shall not pursue a detailed 
calculation of these coefficients for 8 > 1, and shall in- 
stead rely on our previous results for 5 — 1 as a general 



illustration of the potential possibilities for the flavor de- 
pendence. 

These observations are especially relevant for deter- 
mining the neutrino mass eigenvalues resulting from our 
mass matrices. Indeed, one of the unique features of these 
higher-dimensional scenarios involving right-handed neu- 
trinos in the bulk is that it is the same mass matrix 
that determines not only the neutrino flavor mixings but 
also the neutrino mass eigenvalues. This is different from 
the traditional four-dimensional framework in which cou- 
plings to heavy right-handed neutrinos determine the 
neutrino masses (through the traditional seesaw mech- 
anism), whereas mixings amongst the light left-handed 
neutrinos are primarily responsible for determining the 
effective flavor mixing angles. 

It turns out that the eigenvalue phenomenology is very 
different depending on whether S = 1 or 5 > 1. For 5=1, 
we have the mass matrix given in Eq. 12Ufl . However, in 
this case the effective brane/bulk coupling m is necessar- 
ily larger than the effective brane/brane components m a fj 
by a factor (mi?) -1 3> 1. Thus, the mixing with the bulk 
zero mode is dominant and we must actually invert this 
(m + l)-dimensional mass matrix in order to examine its 
mass eigenvalues. We find that in general, Uf — 1 eigenval- 
ues will be of (the smaller) size M 2 a, while the remaining 
two eigenvalues will be of (the larger) size m; likewise, the 
eigenstates corresponding to A ~ M 2 a will have an ex- 
tremely small component involving the right-handed bulk 
neutrino, while the eigenstates corresponding to A ~ m 
will have a fairly significantly component involving the 
right-handed bulk neutrino, with probability w i. Thus, 
the nf — 1 lightest neutrino components are easily inter- 
pretable as our light flavor neutrino mass eigenstates — 
with significant flavor mixings amongst them — provided 
we choose rif to exceed our eventual desired number of 
light neutrinos by one. 

Given these results, we see that our scenario for 5=1 
actually requires an additional left-handed neutrino on 
the brane which participates in the flavor oscillations but 
is otherwise exceedingly heavy. In some sense, this neu- 
trino may be regarded as an additional "sterile" neutrino, 
since its flavor index need not correspond to any of the 
flavors exhibited by the charged fermions of the Standard 
Model; moreover, by choosing an appropriate location 
for this neutrino, its flavor mixings with the other left- 
handed neutrinos can be minimized and/or eliminated. 
However, the unique prediction of this model (as opposed 
to other higher-dimensional models) is that this sterile 
neutrino is necessarily left-handed, living on the brane. 
Of course, its mass is of the same scale (~ m) as that of 
the right-handed zero-mode neutrino living in the bulk. 
Note that the remaining KK bulk excitations have masses 
scaling as ~ nj R, and thus their masses are even heavier 
by an additional hierarchical factor of (mi?) -1 3> 1. 

For 5 > 1, by contrast, the brane/bulk coupling m is 
suppressed by additional factors of (a/R) <§; 1 relative 
to the brane/brane components m a p. Thus, assuming 
that a/R <C (mi?) 2 , we see that the brane/bulk coupling 
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will be extremely small compared to the brane/brane 
components, leading to rif eigenvalues of mass ~ M%a 
and a single eigenvalue of mass ~ <r/R 2 . Although this 
single mass eigenstate will be much lighter than the ri{ 
active eigenstates, it will be overwhelmingly sterile and 
thus will not interact with the Standard-Model particles. 

We see, then, that regardless of the number of extra 
spacetime dimensions, our active neutrino mass eigen- 
states will have masses which scale as ~ M^a in the cur- 
rent setup. There are therefore a number of values for M* 
and a which can yield neutrino masses in a phenomeno- 
logically acceptable range. One appealing possibility, for 
example, is to take M* (the fundamental mass scale in 
our theory) of order M* ~ O(TeV), while our three- 
branes have widths of order a ~ ^(-^pianck)- By con- 
trast, the only solution for which M* and <j~ 1 are of the 
same order of magnitude requires that we take M* in the 
sub-eV range. Likewise, if we demand a ~ ©(TeV^ 1 ), 
then we must take M* ~ C(keV)-C(MeV). However, 
there are other ways of adjusting the overall mass eigen- 
values in this scenario without affecting the flavor mixing 
angles. For example, if there are two extra dimensions 
whose coordinates are y and z, we can easily imagine 
that the Standard-Model brane is located at z = Z\ while 
the right handed neutrino has bulk excitations in the y- 
direction but is localized around z — Z2 . If both localiza- 
tions have common width a, we can retrace our previous 
steps to find that this location "mismatch" in the sec- 
ond extra dimension has the effect of inserting an extra 
suppression factor 



for each entry in the resulting effective mass matrix M. 
While such an overall rescaling of the mass matrix does 
not affect the magnitudes of the flavor mixing angles, it 
directly suppresses the magnitudes of the mass eigenval- 
ues. Thus, in this way, we are free to adjust our mass 
eigenvalues as needed without disturbing the flavor struc- 
ture of our model. Other similar scenarios are also pos- 
sible. 



IV. DISCUSSION 

In this paper, we have presented a "hybrid" model in 
which the split-fermion scenario and the bulk-neutrino 
scenario are joined together. As we have seen, this joining 
has succeeded in yielding a rich spectrum of flavor mixing 
angles without significant fine-tuning, thereby ameliorat- 
ing some of the difficulties inherent in either approach 
alone. 

Needless to say, we have pursued only the most "min- 
imal" approach towards constructing this hybrid model. 
Further possibilities that we have not exploited in- 
clude, for example, possible variable widths of the Gaus- 
sians |17| , the possibility of several right-handed neutri- 
nos, possible orbifold twists [Lq. or possible additional 



flavor-dependent Yukawa couplings of order one. Wc 
have also avoided primordial flavor mixings on the brane, 
and we have taken our brane neutrinos to have vanish- 
ing Majorana masses. While all of these assumptions 
can be relaxed, our main purpose in this paper has been 
to illustrate the range of theoretical possibilities that 
emerge solely from the joining of the split-fermion and 
bulk-neutrino models. 

Likewise, we have not performed a detailed fit of the 
parameters involved, nor have we tried to embed this 
neutrino sector of the Standard Model within a larger 
split-fermion framework that also involves the quarks and 
charged leptons. However, we caution that in order to 
investigate whether such a hybrid model can accommo- 
date all of the observed neutrino-oscillation properties, it 
is generally not sufficient merely to focus on reproducing 
the preferred mass differences and mixing angles often 
quoted in the neutrino literature. This is because such 
neutrino mixing parameters are usually quoted within 
the theoretical context of a simple two- or three-state 
neutrino mixing. In higher-dimensional neutrino mod- 
els, by contrast, one usually faces very complex, multi- 
component neutrino oscillations with a variety of vary- 
ing oscillation lengths and magnitudes (see, e.g., Ref. (3 
for examples of the complexities involved). Indeed, even 
in the small-coupling limit mR <C 1, one typically has 
a four-neutrino (or multi-neutrino) oscillation, depend- 
ing on the number of bulk neutrinos which are present. 
Thus, in order to perform a proper comparison with ex- 
periment, one must go back to the original oscillation 
data, plotted in terms of appropriate baselines, and per- 
form a direct fit between predicted and observed oscilla- 
tion patterns (and suitable time-averages thereof) with- 
out regard for pre-existing theoretical prejudices. This 
too represents a direction for further investigation. 

There are also several important theoretical issues 
which we have not examined in this paper. One press- 
ing issue, for example, concerns how such hybrid settings 
can arise on a fundamental level. Although mechanisms 
of trapping fermions inside a fat brane have been ex- 
amined closely, and although various ways to realize a 
split-fermion scenario seem to be possible, it would be 
interesting to have a more complete picture in which the 
required fermion locations can be explained from some 
deeper principle and/or generated dynamically. 

Nevertheless, we have seen in this paper that joining 
the split-fermion and bulk-neutrino scenarios has pro- 
duced a rich theoretical model which naturally gives rise 
to large neutrino flavor mixing angles without significant 
fine-tuning. Indeed, it is the coupling to the KK modes of 
the bulk neutrino which eliminates the usual fine-tuning 
of the split-fermion scenario, and which produces the rel- 
atively large mixing angles that we have observed. Thus, 
within the context of the split-fermion scenario, we have 
a natural explanation for the somewhat surprising fact 
that quark mixing angles are generally small while lep- 
ton mixing angles are apparently large: while the lepton 
sector contains a right-handed neutrino which is a gauge 
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singlet and can therefore propagate in the bulk, all of 
the fields of the corresponding quark sector carry gauge 
charges and thus must live within the brane. Thus only 
the lepton sector is capable of receiving the enhancement 
to the flavor mixing angles that comes indirectly from the 
couplings to the KK modes of the right-handed neutrino. 
We believe that this observation is critical, and perhaps 
provides one of the strongest motivations for considering 
such higher-dimensional models of neutrino phenomenol- 
ogy. 
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APPENDIX A: PERTURB ATIVE 
DIAGONALIZATION OF THE MASS MATRIX 

The neutrino mass matrix M shows a typical structure 
which is common for seesaw scenarios. We can therefore 
analyze some of its features by similar means. Following 
standard treatments, let us write M. in the block form 



M 



Ml Mr) 



(27) 



where the Ml block corresponds to the brane neutrinos 
as well as the zero mode of the bulk neutrino, where the 
Mr block corresponds to the excited KK modes of the 
bulk neutrino, and where Mjj corresponds to the cou- 
plings between these two groups of states. Assuming 
that these couplings are small (i.e., assuming mR <C 1, 
so that we are in the "perturbative" regime) , we can then 
approximately block-diagonalize this matrix by finding a 
matrix U such that 



U T MU = mL Ji . (28) 

V M R ) V ' 

In general, the solution for U is given by 

V j \j, ^ with k w MdM^ 1 , (29) 



whereupon we find 

Ml ~ M L - MdM^MZ 



M, 



Mr. 



(30) 



Since our original mass matrix is real and symmetric, the 
eigenvalues are necessarily real. 

The next-order corrections have the form 
MjjMft MpMuMft , and an inspection shows that 



these are suppressed by a factor of mR relative to those 
considered above. Likewise, the mixing between the 
excited KK modes and the brane/zero-mode subsector 
is suppressed by a factor of m™ ± i?< 1. 

Note that in the above, we have chosen to group our 
blocks in such a way that the zero mode of the bulk neu- 
trino is joined with the brane neutrinos rather than with 
the bulk-neutrino excited KK states. This is because the 
KK zero mode is the only mode from the KK tower which 
is not heavy, and which therefore fails to decouple. By 
arranging our blocks in this manner, we are therefore "in- 
tegrating out" only the excited KK modes but retaining 
the zero mode in our low-energy reduced mass matrix 

Ml- 

Given the parameters in our specific model, a straight- 
forward calculation gives 



R 

n 



(31) 



whereupon we find 



KMn 



)- 


( m + 


°) 


1 n 


V rn_ 





-m a p 




(32) 



with 



n=l 



We thus find 



Mr 



R p 



R 



m 




(33) 



(34) 



whereupon inserting the definitions in Eq. Q yields the 
result given in Eq. (jT7|l . 



APPENDIX B: MASS SCALES FOR 8 > 1 

We now consider the mass scales involved in our mass 
matrix M. We are particularly interested in the mass 
scales associated with the expression m a p in Eq. 133|l . 
Thus, there are three components of this expression that 
we need to generalize to 5 dimensions: we need to gen- 
eralize the couplings m" ± ; we need to generalize expres- 
sions such as 1/n; and we need to generalize the KK sum. 
We shall therefore begin by addressing each of these in 
turn. 

For excitations of the bulk field in more than one ex- 
tra spacetime dimension, the Kaluza-Klein expansion in 
Eq. J2Jl will now involve products of sines and cosines. 
This in turn means that the couplings m" ± in Eq. Q 
will also be products of sines and cosines whose argu- 
ments display the positions yi of our fields. 
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Likewise, the form of the higher-dimensional Clifford 
algebra implies that a Dirac spinor in 4 + 6 spacetime 
dimensions will have 2 2+5 / 2 components, and the inverse 
of the wave vector n will be a linear combination of its 
components. However, since our goal is merely to esti- 
mate relevant mass scales, we shall ignore these spinorial 
details and retain a notation whereby we treat quantities 
such as 1/n as a single number. 

Finally, we observe that the KK sum in Eq. i|33|) 
will now extend over all ^-dimensional vectors n = 
(n 1 ,n 2; ■■■ins) with rii 6 2 + except the zero vector. 

Let us now put these ingredients together. First, we 
notice that in the higher-dimensional case with 5 > 1, 
the couplings from Eq. ® will enter into Eq. Ijl5|l only 
for mixtures of even and odd fields. This means that the 
factors that must be summed as part of the n-summation 
will consist of at least one sine function and one cosine 
function. Since the contributions to this sum are natu- 
rally damped beyond n ~ R/cr, and since we assume the 
locations y a to be of order a, we can approximate each 
occurrence of the sine function by the linear term in its 
Taylor expansion (and set each cosine factor to 1). Let 
s denote the number of sine functions in the product, 
where 1 < s < S, and likewise let i = 1, s denote those 
directions for which these sine functions appear. Each 
of these sine functions will then make a linear contribu- 



tion to m a f3 of the form yin/R, whereupon we find that 
Eqs. I|13fl and (|15fl will generalize to yield an expression 
of the form 

m aP ~ m 2 Rj2~^(-^r) f[(w5) ■ ( 35 ) 

n x 7 i— 1 

Note that this n-summation continues to represent a sum 
over ^-dimensional vectors n. Defining k = (a / R)n, we 
see that this sum can be approximated by an integral 

m ap ~ ^JJ^/ff) J dkk s+s ^exp(-k 2 /n 2 ) . 

(36) 

However, the integral is now purely dimensionless and ~ 
0(1); likewise, for j/j ~ a, we see that each of the product 
factors are also ~ 0(1). Thus the overall mass scale is set 
by the prefactor in Eq. (|36[) . However, given the result in 
Eq. 124|) . we see that we can write this prefactor in the 
form M 2 a, obtaining the same overall mass scale as we 
found for 8 = 1. Thus, we conclude that the overall mass 
scale for m a p is independent of both R and <5. Indeed, it is 
easy to verify that this result remains true even when all 
of the proper numerical and spinorial factors are included 
in the analysis. 
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